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Abstract: Models of period variations are basic tools for period ana- 
lyzes of variable stars. We introduce phase function and instant period 
and formulate basic relations and equations among them. Some simple 
period models are also presented. 
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1. Motivation 

The majority of variable stars change brightness (also radial velocity, mag- 
netic field, etc.) more or less regularly and with one period P. As a rule the 
form and amplitude of their light curves remain constant for plenty of cycles 
while their periods may change slightly for a number of reasons. The goal 
of this paper is to create models of such period variations which can then 
serve as basic tools for the advanced period analysis of changes of variable 
objects. 



2. The instantaneous period and the phase function 

The state of a periodically changing variable star is described by two functi- 
ons of time t: the non-descending stairs-like epoch function E(t) expressing 
the number of cycles elapsed from the moment beginning of epochs counting 
Mq, and the sawtooth function called the phase '^{t), reaching its minimum 
{'^{t) = 0) when the new cycle of the variability starts and its maximum 



((/^(t) = 1) when the particular cycle ends. The phase is used for the con- 
struction of so called phase curves of various values characterizing the va- 
riable star. The instantaneous period for a selected epoch P{E) is then 
the duration of this epoch. It is useful to introduce a new time-dependent 
quantity - phase function as 

^{t) = E{t) + if{t); ip{t) = frac[i?(t)]; E{t) = floor[i9(t)], (1) 

that could substitute both these 'mathematically awful' functions. The ope- 
rator 'frac' removes from a real number its integral part, while 'floor' rounds 
the quantity to the nearest lower integer. 

The phase function 'd{t) is a monotonic rising function of time originating 
at the beginning of epoch counting t = Mq, 't?(Mo) = 0. Using we are 
able to determine the epoch and the phase (see Eq.[T]) for any time. The time 
derivative of the phase function at the time t is equal to the instantaneous 
frequency i}{t) = f{t) = P{t)~^, where P{t) is the instantaneous period of 
the variable star at the time t (firstly introduced in [1]). The phase function 
is then determined by the following differential equation 

""W- 1 . ,(, = M„)=0; ^ m=t 4^- (2) 



it P(ty ^ • " iM„P(r) 

It is useful to also introduce the inversion function T(i?) to the phase 
function {} which determines the moment T when the phase function rea- 
ches the assigned value. The times of the zero phase (usually the times of 
brightness minima or maxima) O for the distinct epoch E are then given 
by the relation: e{E) = T{i} = E). 

The function T(??) can be derived as the inversion function of ■i?(t) or it 
can be found by solving the differential equation: 

dr(??) r-9 



di? 



P{{}); T{S = 0) = Mo; ^ T{S) = + f P{C) dC (3) 

Jo 



3. Some simple period models 

Let us derive the phase functions, their inverses and period time dependen- 
cies for some astrophysically important polynomial period models. 

P{t) = Pq 

The simplest period model that we use as a first approximation assumes that 
the period of variations is constant: P{t) = Pq. Then, using the equations 
([2]) and da]), we get for its t?i(t) and Tiji?) 

Mt) = *-^; Ti(t?) = Mo + Po^i. (4) 
Pq 



P{t) = PQ 

Let's assume that the period is a hnear function of the time: 



P = Po; P{t) = Po + Po{t-Mo) = Po{l + Po^i), (5) 
dP_dP_dPd^9_dPl _ . . p.., p P„,. ... 



^(t) = ^ln(l + Po^i); T(^) = Mo + S (e^o^ - l) . (7) 
Pq Po 



Since the instantaneous rate of period P is as a rule very slow, we can 
replace the real phase functions ■d{t) and their inverses by their Maclaurin 
decomposition: 

m = ^ln(l + Po^i) = ^i-Pof + Po't 



T(^) = Mo + -^(e^»^-lj=Mo + Po(^^ + Po^ + Po'^--j- (9) 

For the majority of variable stars with an inconstant period it is valid that 
P At <C P, where At is the whole duration of its observation. The last 
decomposition terms in (see Eqs. El E]) can then be neglected. 

= - f (^)^ rW . M„ + P.. + PoP^; 

= Po + ^'o (t - Mo); P(i?) = Po (l + Po • (10) 

The phase function and its inverse are quadratic functions and to describe 
then we just need three parameters: Mo, Pq, Po- 

The exact same relations as given in Eqs. (|10|) can be obtained if we 
assume other laws for period development, such as PP = PoPo or P/P = 
Pq/Pq, or if we assume the parabolic course of the phase function i?(t) or 
its inverse ^(^?), T{E). 

P = Pn 



If we find that the dependence of i9(t) have apparent cubic or higher terms, 
we acknowledge that the first derivative of the period is inconstant. Then 
the simplest period model assumes the period to be a quadratic function oi 
time: 



,2o ^? dPd^9i 



P{t) = Po + Po Po + P^Po y; P{t) = ^^^ = Po + PoPo ^1. (11) 



Using Eq. ([2]) we can find the exact solution for the phase function 'd{t). 
However, the function expression is so complex that it is practically useless. 
Consequently we only the first three terms of its Mclaurin decomposition. 



„q2 „q3 ,q2 q3 

^(i) = ^i-^o^ -(^'o^o-2Po')^ = ^i-^o^ -PoPo^; (12) 

„q2 „q3 
r(i?) = Mo + Po^ + PoPo^ + {Pi Po + Po P^)^ 

= Mo + Po^ + PoPo^ + PiPo^; (13) 
= 1 + ^ + (Po Po + p2) = 1 + Po ^ + Po Po (14) 



The expressions in brackets in Eqs. p2|) . p3|) . and (|M|) can be simplified, 
because of PqPo S> Pq is valid for all known cases of variable stars changing 
their period. 



4. Conslusion 

We introduced several simple polynomial period models that can generally 
be used for most known and unknown period trends. More complicated 
models of period changes in eclipsing binary systems, accounting for the light 
time effect caused by third bodies as well as apsidal motion, should also be 
developed using the relations given in Section 2. However, the formulation 
of those period models is beyond the scope of this article. 
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